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Abstract 

In this paper the hp-version of the boundary element method is apphed to the electric 
field integral equation on a piecewise plane (open or closed) Lipschitz surface. The under- 
lying meshes are supposed to be quasi-uniform. We use H(div)-conforming discretisations 
with quadrilateral elements of Raviart-Thomas type and establish quasi-optimal convergence 
of /ip-approximations. Main ingredient of our analysis is a new H^^/^(div)-conforming p- 
interpolation operator that assumes only H'' n H~"'^/^(div)-regularity (r > 0) and for which 
we show quasi-stability with respect to polynomial degrees. 
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1 Introduction and formulation of the problem 

In this paper we prove convergence of the hp-version of the boundary element method (BEM) for 
the electric field integral equation (EFIE) on piecewise plane (open or closed) surfaces discretised 
by quasi-uniform meshes. The EFIE is a boundary integral equation that represents a boundary 
value problem for the time-harmonic Maxwell's equations in the exterior domain. It models the 
scattering of electro-magnetic waves at a perfectly conducting body (the scatterer). The solution 
to the EFIE is the induced electric surface current (a tangential vector field) on the surface of the 
scatterer, see, e.g., [28j . If the scatterer is a thin object (i.e., its thickness is small in comparison 
to the wave length), then it can be modelled as an open surface (a sub- manifold with boundary) 
in M^. Our analysis covers this theoretically challenging case, which has important applications 
(e.g., antenna problems). 

The basis of our BEM is a variational formulation of the EFIE, called Rumsey's principle. 
For smooth surfaces, its boundary element discretisation has been studied by Bendali in 1984, see 
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[21 S]. Progress in the numerical analysis of the EFIE on Lipschitz surfaces has been achieved 
relatively recently and was inspired by the study of traces of functional spaces that govern 
Maxwell's equations in Lipschitz domains [1^. The main challenges in this analysis concern 
the solvability and quasi-optimal convergence of approximations and a priori error estimation 
in the energy norm. In the framework of the /i- version of the BEM, i.e., for discretisations 
with elements of fixed order on refined meshes, these issues were addressed in [15^ [2^ [T8] (for 
polyhedral surfaces) and in [12] (for open Lipschitz surfaces). We note that in [231 UHl [l2] 
the authors focused on conforming discretisations of Rumsey's principle with Raviart-Thomas 
(or Brezzi-Douglas-Marini) boundary elements, whereas in |15j a mixed formulation utilising 
standard (continuous) basis functions was used. In this paper we follow the former approach, 
called the natural boundary element method for the EFIE. 

While in the p-version of the BEM the mesh is fixed and approximations are improved by 
increasing polynomial degrees, the hp-version combines both mesh refinement and increase of 
polynomial degrees. In our previous paper ^ we analysed the natural p-BEM for the EFIE on a 
plane open surface with polygonal boundary. We have proved convergence of the p-version with 
Raviart-Thomas (RT) parallelogram elements and derived an a priori error estimate which takes 
into account the strong singular behaviour of the solution at edges and corners of the surface. 
With the present note we prove the convergence of the natural /ip-BEM on polyhedral and 
piecewise plane open surfaces discretised by quasi-uniform meshes of quadrilateral (in general, 
curvilinear) elements. We emphasize that RT-spaces are used for both affine and non-affine 
quadrilateral meshes. 

In order to prove convergence of approximations for the EFIE, one usually relies on properties 
of the continuous and discrete Helmholtz decompositions, and on the proximity in some sense 
of the discrete decompositions to the continuous one, see [I2l HI]. The key property is the or- 
thogonality of decompositions, and the main tool in the analysis is an appropriate interpolation 
operator (a projector onto the corresponding polynomial space). In [121 [18], L^-orthogonal dis- 
crete decompositions mimicking the Helmholtz decomposition of the energy space were analysed 
for finite dimensional subspaces based on Raviart-Thomas and Brezzi-Douglas-Marini (BDM) 
boundary elements. It has been proved that these discrete decompositions are sufficiently close 
to the continuous one as the mesh parameter h tends to zero (i.e., for the /i- version of the BEM). 
The main tools in the proofs were the standard H(div)-conforming RT or BDM interpolation 
operators. However, it turns out that a generalisation of that approach to the p- and the hp- 
versions is not straightforward when sticking to both the L'^-orthogonality of decompositions and 
classical interpolation operators. In [8j, for the p- version, we employed an H~^/^-orthogonality 
of the Helmholtz decomposition, while using the classical RT interpolation operator. Unfortu- 
nately, the extension of this approach to polyhedral surfaces and even to piecewise plane screens 
does not seem to be easy. In particular, the low regularity of the Laplace-Beltrami operator 
on polyhedral surfaces is not enough to prove stability (with respect to polynomial degrees) 
of the classical Raviart-Thomas interpolation operator. That is why in this paper we use an 
alternative approach: we adhere to the L^-orthogonality of decompositions but utilise a non- 
classical H~^/^(div)-conforming interpolation operator (for RT-elements). Our construction of 
this operator (see Section [5|) is much in the spirit of [25 > where H^- and H(curl)-conforming 
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projection-based interpolation operators were introduced and analysed. We, however, employ 
the ^^^/^-projection for the divergence term, and thus need the corresponding inner product to 
be written in an appropriate explicit form. Moreover, the use of an appropriate scaling argument 
allows to prove the convergence result in the framework of the /ip-BEM. 

We will denote by F a piecewise plane (open or closed) Lipschitz surface in M^. Let us 
introduce Rumsey's formulation of the electric field integral equation on T. For a given wave 
number k > and a scalar function v (resp., tangential vector field v) we define the single layer 
operator ^'^ (resp., by 

^kvix) = — v{y)j r dSy, X G M^F 

4vr \x -y\ 

fresp., *fcv(x) = ^ / v(y) ^ ' ^' dSy, x G M^f). 

Let L^(F) be the space of two-dimensional, tangential, square integrable vector fields on F. 
By Vr (resp., divr) we denote the surface gradient (resp., surface divergence) acting on scalar 
functions (resp., tangential vector fields) on F. We will need the following space: 

X = H-V2(divr,F) := {u G H^^/^(F); divr u G R~^I'^[V)\ 
if F is a closed surface, and 

X = Ho^/^(divr,F) := {u G H[^^^(F); divr u G ^-^/^(F) and 

(u, Vru) + (divr u, i;) = for all v G C°°(F)} 

if F is an open surface. In the latter definition the brackets (•, •) denote dualities associated with 
H|'|''^(F) and //^/^(F), respectively. For definitions of the space C°°(F) and the Sobolev spaces 
on F see ^3.11 

Throughout, we use boldface symbols for vector fields. The spaces (or sets) of vector fields 
are also denoted in boldface (e.g., H^(F) = (i7*(F))'^), with their norms and inner products 
being defined in §3.11 

Let X' be the dual space of X (with (F) as pivot space). Now, for a given tangential vector 
field f G X' (f represents the excitation by an incident wave), Rumsey's formulation reads as: 
find a complex tangential field u G X such that 

a(u,v) := (7to(^fedivru),divrv) -A:2(,r,(*^u),v) = (f,v) Vv G X. (1.1) 

Here 7tr is the standard trace operator, and vr^ denotes the tangential components trace mapping 
(see §3.11 for definitions and properties of these operators). 

The rest of the paper is organised as follows. In the next section we define the /ip- version of 
the BEM for the EFIE and formulate the main result (Theorem 12. ID . which states the unique 
solvability and quasi-optimal convergence of this approximation method. Section [3] gives nec- 
essary preliminaries: in §3.11 we recall definitions of functional spaces of scalar functions and 
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vector fields; then in ^3.21 we introduce some equivalent norms in the Sobolev spaces and 
H'^ {s = ib^) on the reference element and derive expressions for corresponding inner products; 
some auxiliary lemmas are collected in ^3.3i In Section S] we discuss the continuous Helmholtz 
decomposition of the energy space X and its discrete counterpart. Section [5] is devoted to inter- 
polation operators. First, we recall some known operators such as the L?- and ^^^/^-projectors, 
the H^- and H(curl)-conforming projection-based interpolation operators. Then we introduce 
an H~^/^(div)-conforming interpolation operator and study its properties. In particular, we 
prove the quasi-stability of this operator and its commutativity with the -ff~^/^-projector. The 
H~^/^(div)-conforming interpolation operator and L^-orthogonal discrete Helmholtz decompo- 
sitions are the main tools in the proof of Theorem 12. H which is given in Section [6l We conclude 
the paper with Section [71 where we comment on extensions and open problems. 

Throughout the paper, C denotes a generic positive constant which is independent of /i, p, 
and involved functions, unless stated otherwise. 

2 The version of the BEM and the main result 

For the approximate solution of (jl.ip we apply the /ip-version of the BEM based on Galerkin 
discretisations with Raviart-Thomas spaces on quasi-uniform meshes. In what follows, /i > 
and p > 1 will always specify the mesh parameter and a polynomial degree, respectively. For 
any Q. C we will denote po, = sup{diam(i?); B is a ball in 0}. Furthermore, we will denote 
hy K = (0, 1)^ the reference square. The sides of K will be denoted by (i = 1, . . . , 4). 

Let T = {A/i} be a family of meshes A/^ = {Fj; j = 1, . . . , J} on F, where the elements Fj are 
open quadrilaterals (in general, curvilinear ones) satisfying the following standard assumptions 
(below we denote hj = diam(Fj) for any Fj G A/j): 

i) F = U^^j^f J-; the intersection of any two quadrilaterals Fj, F^ {j ^ k) is either a common 
vertex, an entire side, or empty; 

ii) The elements are shape regular, i.e., there exists a positive constant C independent of 
h = max hj such that for any Fj E Ah and arbitrary Ah G T there holds hj < C pr^ ■ 

iii) Any element Fj is the image of the reference square K, more precisely 

f j = Tj{K), x = rj(0, X = {xi,X2) G f,, ^ = (6,6) e K, 

where Tj is sufficiently smooth one-to-one mapping with sufficiently smooth inverse T^^^ : 
Fj K. The Jacobian matrix of Tj is denoted by DTj{^), it is supposed to be invertible 
for any K, and DTf^i^x.) = {DTj{$))-'^. We assume that 

|det(Drj(^))| ~ h] for any ^£K, (2.1) 
and there exist positive constants C independent of h such that for k = 1,2 there holds 

sup \\D%m\c,[w^m < ch';, sup ||z)'=r-i(x)ll£,(R2,M2) < chj\ (2.2) 
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Here D^Tj(^) (resp., D^Tj ^{'x.)) denotes the k-tli (Frechet) derivative of Tj (resp., Tj ^) 
and II • ||£j.(x,y) is the operator norm in the space Ck{X, Y) of continuous fc-hnear mappings 
from into Y . 

iv) if Tj n f is an entire side then denoting by (resp., r|) the restriction of Tj (resp., 
Tfe) to the corresponding side T~^{£) (resp., r^^(^)) of the reference element K, one has 
Tj = Tf. as mappings of the unit interval (0, 1) onto i. 

Remark 2.1 Assumptions (|2.ip . (|2.2p above are always satisfied for affine families of elements 
(i.e., for meshes of parallelograms). For curvilinear elements these assumptions are satisfied, 
for instance, if the elements tend to he affine as h ^ (see \19\ Section 4.3]j. In this case ()2.ip . 
(|2.2p follow from assumption ii) and the smoothness ofTj, provided that hj is small enough. 

In this paper we consider a family T of quasi-uniform meshes Ah on F in the sense that 
there exists a positive constant C independent of h such that for any Tj G Ah and arbitrary 
Ah € T there holds h < C hj. 

The mapping Tj introduced above is used to associate the scalar function u defined on the 
real element Tj with the function u defined on the reference element K: 

u = uo T^^ on Fj and u = uoTj on K. 

Any vector-valued function v defined on K is transformed to the function v on Fj by using the 
Piola transformation: 

V = Mj{^) = j-DTjV o Tr\ v = Mj\v) = JjDTr^v o Tj, (2.3) 

where Jj is the determinant Jj{^) ■= det{DTj{^)). 

Let us introduce the needed polynomial sets. By Vp{I) we denote the set of polynomials 
of degree < p on an interval / C M, and Vp{I) denotes the subset of Vp{I) which consists of 
polynomials vanishing at the end points of /. In particular, these two sets will be used for edges 
£i C dK. 

Further, 'Pp^^p^{K) denotes the set of polynomials on K of degree < pi in and degree 
< p2 in ,^2- For pi = p2 = p we denote Vp{K) = Vp^p{K). The corresponding set of polynomial 
(scalar) bubble functions on K is denoted by Vp{K). 

Let us denote by V^iK) the RT -space of order p > 1 on (see, e.g., [29]). i.e., 

Vj'iK) =Vp,p-i{K) yiVp-i,p{K). 

The subset of {K) which consists of vector- valued polynomials with vanishing normal trace 
on the boundary dK (vector bubble-functions) will be denoted by ^^^'^{K). 
Then using transformations (j2.3p . we set 

^hp := {v G X"; A^J^vlrJ G Vf^{K), j = 1, . . . , J}, (2.4) 
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where the space X*^ C X is defined in ^3.11 (X*^ = H(divr, T) if F is closed and X'^ = Ho(divr, F) 
if F is an open surface). We will denote by = N{h,p) the dimension of the discrete space 
X/ip. One has ~ h~'^ for fixed p and N c^p^ for fixed h. 

The /ip-version of the Galerkin BEM for the EFIE reads as: Find u^p G X/^p such that 

aKp,v) = (f,v) VvGX;,p. (2.5) 

Let us formulate the result which states the unique solvability of (|2.5p and quasi-optimal 
convergence of the /ip- version of the BEM for the EFIE. 

Theorem 2.1 There exists Nq > 1 such that for any f G X' and for arbitrary mesh- degree 
combination satisfying N{h,p) > Nq the discrete problem (12. Sp is uniquely solvable and the 
hp-version of the Galerkin BEM generated by RT-elements converges quasi- optimally, i.e., 

||u - u,,p||x < Cinf{||u - v||x; v G X/,p}. 

Here, u G X zs the solution of (jl.ip . u/jp G X/jp is the solution of ()2.5p . || • ||x denotes the norm 
in X, and C > is a constant independent of h and p. 

The proof of Theorem 12.11 is given in Section [6] below. 



3 Preliminaries 

3.1 Functional spaces, norms, and inner products 

First, let us recall the Sobolev spaces and norms for scalar functions on a Lipschitz domain 
Q C M", see |25) . To that end we will need the space C°°{il.) of infinitely differentiable functions 
in S7 and its subspace CQ°(r2) C C°°(rj) which consists of functions with compact support in fi. 
For an integer s, let H'^{Q) be the closure of C°°(0) with respect to the norm 

where 



[ \D'u{x)fdx, and H^{n) = L'^{n). 
Jn 



Here, \D^u{x)\'^ = X]|o|=s \D°^u{x)\'^ in the usual notation with multi-index a = (oi, . . . ,««) and 
with respect to Cartesian coordinates x = (xi, . . . ,Xn)- For a positive non-integer s = m -\- a 
with integer m > and < a < 1, the norm in H^{il.) is 

II l|2 — II ||2 _|_ I l2 

with semi-norm 

\a\=m 
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The Sobolev spaces H^{0,) for s S (0, 1) and for a bounded Lipschitz domain are defined by 
interpolation. We use the real K-method of interpolation (see |25j) to define 

H'{n) = (^L'^{Q),H*q{Q)^ ^ ^ (1/2 < t < 1, 0<s<t). 

Here, F^(17) (0 < t < 1) is the completion of C^{Q) in H\Q) and we identify i?o(^) and 
H^{n). Note that the Sobolev spaces H'^{0,) also satisfy the interpolation property 

H%n) = (I'^in), H^ifi)) (0 < s < 1) 

with equivalent norms. Furthermore, the semi- norm | • {^^-{n) is a norm in H^{n) due to the 
Poincare inequality. For the L^(0)-norm we will use the notation || • ||o,r2- 

For s G [—1,0) the Sobolev spaces and their norms are defined by duality with L^(ri) = 
H^{n) = H^{n) as pivot space: 

\MH-in)= sup -—r , 1^11^.(0)= sup — , (3.1) 

Oytv£H-={n) W^Wn-^iu) OjtveH--{n) \\v\\H--{n) 



where 



{u,v) = {u,v)o^Q := / u{x)v{x)dx 
Jn 



denotes the extension of the L^(J7)-inner product by duality (and v is the complex conjugate of 
v). 

Now, let r be a piecewise smooth (open or closed) Lipschitz surface in M^. We will assume 
that r has plane faces F^*) {i = 1, ... ,1; without loss of generality it is assumed that T > 1) 
and straight edges etj = is a closed surface, we will denote by 17 

the Lipschitz polyhedron bounded by F, i.e., F = dCl. In the case of an open surface F, we first 
introduce a piecewise plane closed Lipschitz surface F which contains F, and then denote by 
Q the Lipschitz polyhedron bounded by f , i.e., f = dil.. For each face F^*) C F there exists a 
constant unit normal vector i/j, which is an outer normal vector to 0. These vectors are then 
blended into a unit normal vector u defined almost everywhere on F. For each pair of indices 
i,j = I, ... ,2 such that f n f = Cij we consider unit vectors Tij, Tp\ and r^*^ such that 

Tij\\eij, r'f^ = Tij X Ui, and r^*^ = Tij x Uj. Since each F^*) can be identified with a bounded 

subset in M^, the pair (rp^Tjj) is an orthonormal basis of the plane generated by F^*^ 

Let F be a closed surface. Then F = is locally the graph of a Lipschitz function. 
Since the Sobolev spaces H'^ for \s\ < 1 are invariant under Lipschitz (i.e., C^'^) coordinate 
transformations, the spaces H^{T) with |s| < 1 are defined in the usual way via a partition of 
unity subordinate to a finite family of local coordinate patches (see [26l)- Due to this definition, 
the properties of Sobolev spaces on Lipschitz domains in carry over to Sobolev spaces on 
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Lipschitz surfaces. If T is an open surface, then the Sobolev spaces H'^(T), H^(T) for |s| < 1 
and Hq{T) for < s < 1 are constructed in terms of the Sobolev spaces H^{T) on a closed 
Lipschitz surface f D T (see [26]). Note that the spaces iJ*(r(*)) and H^{r^^'^) on each face T^*) 
are well-defined for any s > — 1. 

For s > 1 we define the space H^(T) in the following piecewise fashion (hereafter, Ui denotes 
the restriction of u to the face F^*^): 

H%T) := {u E H\T); G i/^(F«), i = l,...,I}. 

This space is equipped with its natural norm 



J 

l2 II ||2 



\H^r) ■— I ll«'ll//i(r) "T ll"'«llii-'i(r(')) 
^ i=l 

Besides the above, we will need the following spaces: 

H:{T) :={ueH'{Ty, (u,l)o,r = 0}, 

where s > — 1 if F is closed, and s > — ^ if F is an open surface. 

We will denote by 7tr the standard trace operator, jtriu) = u\r, u £ C°°(fi). For s G (0, 1) 
(resp., s > 1), 7tr has a unique extension to a continuous operator i?*"'"^/^(r2) — > H^{T) (resp., 
H'+^/^{n) H^{T)), see [SnillS]. We wiU use the notation C°^{r) = 7tr(C°°(fl)). 

Using the introduced Sobolev spaces of scalar functions, we define for s > — 1: 

H"(0) = iH'{n)f, H"(F) = iH'{T)f; 

H"(rW) = (i7«(FW))2, ff(F(^)) = (^^(fW))2, 1 < i < Z. 
If F is an open surface, then in addition to the above we define the space 

H^(F) = {H'{T)f, \s\ < 1. 

The norms and inner products in all these spaces are defined component- wise and usual conven- 
tions HO(rj) = L2(17), H0(F) = H0(F) = L2(F), H0(F«) = H0(F«) = L2(F«) hold. 

Now let us introduce the Sobolev spaces of tangential vector fields defined on F (see [131 El 
[16]). We start with the space 

L^{r) := {u G L2(F); u • = on F}, 

which will be identified with the space of two-dimensional, tangential, square integrable vector 
fields. The norm and inner product in this space will be denoted by || ■ ||o,r and (•, •) = (•, •)o,r) 
respectively. The similarity of this notation with the one for scalar functions should not lead to 
any confusion. Then we define: 

H^_(F) := {u G L2(F); Ui G H^(F«), 1 < i < T}, s > 0, 
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l|u||H''_(r) :- ( E lluilln'^frW) 

\ i=l 

Let 7tr be the trace operator (now acting on vector fields), 7tr(u) = u|r, 7tr : H'^^^/^(r2) — > 
H*(r) for s G (0, 1), and let 7^^^^ be one of its right inverses. We will use the "tangential 
components trace" mapping tTt- : {C°°{^))^ L^(r) and the "tangential trace" mapping 7,- : 
(C°°(ll))^ — > Lj(r), which are defined as u 1-^ x (u x i/)|p and u x i/|p, respectively. We will 
also use the notation iTr (resp., 7,-) for the composite operator tt^ o 7"^ (resp., 7,- o 7^^"*^), which 
acts on traces. Then we define the spaces 

Hj/^r) := 7r.(HV2(r)), nf (F) := jAU'/^T)), 

endowed with their operator norms 

ll^lln^/^rn •= ^ ™! {ll</'llHV2(r); 7rr(0) = u}, 
"11 (?!)eHi/2(r) ' 

l|u||„i/2.ps := inf {||0||Hi/2(r); 7r(0) = u}. 

1/2 1/2 

It has been shown in [13] that the space Hy (F) (resp., Hj^ (F)) can be characterised as 

1/2 

the space of tangential vector fields belonging to H _ (F) and satisfying an appropriate "weak 
continuity" condition for the tangential (resp., normal) component across each edge Cij of F. 
For s > I we set 

H||(F) := {u G H*_(F); Uj • tij = uj ■ Tij at each Cij}, 

Hi(F) := {u G H^_(F); u, • r^^ = u^- • rf at each dj}. 

For any s > ^ the spaces Hjj(F) and H5_(F) are closed subspaces of H*_(F). Finally, for s G [0, ^) 
we set 

H^(F) = Hi(F) := H^_(F). 

If F is an open surface, then we also need to define subspaces of H|(F) and H^(F) incorpo- 
rating boundary conditions on dT (for tangential and normal components, respectively). In this 
case, for a given function u on F, we will denote by u the extension of u by zero onto a closed 
Lipschitz polyhedral surface F D F. Then we define the spaces 

Hjj(F) := {u G H|(F); a G H|(f)}, s > 0, 

Hi(F) := {uGHi(F); liGHi(f)}, s>0, 
which are furnished with the norms 

ll"llH^(r) ll"llHj|(f)' ll"llH5^(r) ll"llH^(f)' s>0. 
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When considering open and closed surfaces at the same time we use the notation Hj|(r), H^(r), 
etc. also for closed surfaces by assuming that Hjj(r) = Hjj(r), H^(r) = H^(r), etc. in this 
case. This in particular applies to the following definition of dual spaces. For s G [—1,0), the 
spaces H|(r), H|(r), H^(r), H^(r), and H"_(r) are defined as the dual spaces of H-"(r), 

H-^(r), H;^(r), Hj'(r), and H-'(r), respectively (with L'f{T) as pivot space). They are 
equipped with their natural (dual) norms. Moreover, for any s € ( — ^) there holds (cf. [22]) 

H|(r) = H|(r) = Hi(r) = Hi(r) = ff_(r). 

Using the above spaces of tangential vector fields, one can define basic differential operators 
on r. The tangential gradient, Vr : H^(r) L^(r), and the tangential vector curl, curlp : 
H^{T) L?(r), are defined in the usual way by localisation to each face T^^^ (see [Ml E] for 
definitions and properties of these operators on both closed and open surfaces). To proceed with 
extensions of these operators and with their adjoints we need to distinguish between open and 
closed surfaces. 

Let r be a closed surface. The adjoint operators of — Vr and curlp are the surface divergence 
and the surface scalar curl, respectively: 

divr:L2(r)^/7-i(r), cm\r : Lf{T) ^ H;\T). (3.2) 

It has been shown in [14J that Vr and curlr can be extended to 

Vr : H^/^{r) ^ H;^/^(r), curlr : H^^^{r) ^ H~^/^(r). 
Moreover, they have closed ranges in corresponding spaces. Their adjoint operators 

divr : Hy'(r) ^ H;'/\t), curlr : Hj/'(r) ^ H;'^\t) 

are linear continuous and surjective. 

Finally, the Laplace-Beltrami operator is defined on T as follows 

Aru = divr(Vr u) = - curlr (curlr u) Vn e -H'^(r). (3.3) 

One has Ar : H^(r) — > H^^{T), it is linear continuous and invertible. 
If r is an open surface, then instead of ()3.2p there holds 

divr : L^{r) i?-^(r), curlr : L2(r) ^ -H'"^(r). 

The operators Vr and curlr again can be extended as follows (cf. [14J): 

Vr : H^^^iT) ^ H;^/^(r), Vr : H^/^{r) ^ H7^/^(r) 

and 

curlr : H^^Hr) ^ H"^/^(r), curlr : H^^^{r) ^ H"^/^(r); 
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they also have closed ranges in corresponding spaces, and their adjoints 

divr : uY\r) ^ H-^/\T), divr : Hy^(r) ^ H~^/\T) 

and 

curlr : Hj/^(r) ^ R-^/^T), curlr : Hj^^(r) ^ H~^^\T) 

are linear continuous and surjective. The Laplace-Beltrami operator Ap : H^(r) H^^{T) is 
defined as in (j3.3p . 

We will need the following spaces involving Ap: 

7^(r) := {u£H^{r)/C; ArueH~^^^{r)} if T is closed, 
n{T) ■= {u(iH^{T)/C; ATueH-^''^{T) and 

(Vr u, Vr v) + (Ap u,v) = Q for all v G //""^(F)} if F is an open surface. 

Now we can introduce the spaces which appear when dealing with the EFIE on F. First, we 

set 

H"(divr,F) := {u G H^(F); divp u G H'{T)}, s G [-1/2,0] 

(here, F is either a closed or an open surface). If F is an open surface, then we will also use the 
space 

H^(divr,F) := {u G H|(F); divr u G H'{T)}, s G [-1/2,0]. 

The spaces H*'(divr,F) and H*(divr,F) are equipped with their graph norms |[ • ||H''(divr,r) 

II • llHs(divrr)' respectively. For s = we drop the superscript and for open surfaces also the 

tilde in the above notation, H°(divr,F) = H°(divr,F) = H(divr,F). 

On open surfaces, one also needs the spaces incorporating homogeneous boundary conditions 

~ —1/2 

for the trace of the normal component on dT. By Ho(divr, F) (resp., Hg (divr, F)) we denote 
the subspace of elements u G H(divr, F) (resp., u G H~^/^(divr, F)) such that for all v G C°°(F) 
there holds 

(u, Vrv) + (divr u, v) = 0. (3.4) 

We note that if u G Hq ^^^(divr, F), then identity (|3.4p holds for any v G H^/'^{T) by density. In 
particular, Hg ^''^(divr,F) is a closed subspace of H~^/'^(divr, F). To join the notation for open 
and closed surfaces, we will write 

X° = H(divr,F), X = 11-^/^ {divr, T) if F is a closed surface, 

X'^ = Ho (divr, F), X = Hq ^''^(divr, F) if F is an open surface. 
In both cases the norm in the space X will be denoted as || • ||x- 
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3.2 Some equivalent norms and corresponding inner products 

In this subsection we consider the Sobolev spaces and on the reference square K for 
s = ib^. We will derive expressions for norms which are equivalent to those introduced in 
Section 13. 1[ We note that all results of this subsection are valid also for K being the equilateral 
reference triangle. First, let us introduce some notation. 

1°. We denote by Vt the cube Vt = K x (0,1). Thus = uf^^fi. Let K = Ti _= 
{(2;i,X2,0); {xi,X2) G K}, = {(xi,rE2,l); (xi,X2) G K}, and denote K = dQ}^^. 
Note that K is an open surface. We will use the standard notation for the gradient V and 
for the Laplace operator A, both acting on scalar functions of three variables. 

2°. Given u e H^^^'^{K), we denote by uk the solution of the mixed problem: find uk G 
H'^{n) such that 

Auk = Omn, ^ = uonK, ur = on dn\K. 

If ti G H^^^^{K), then we will use the same notation as above with K replaced by K. 

3°. Given u G H^^'^{d^l), we denote by u its harmonic extension, i.e., the solution of the 
Dirichlet problem: find u G //^(O) such that 

An = in 17, u = u on d^l. (3-5) 
4°. Given u G H^^'^{K), we denote by u° the extension of u by zero onto dO,. Thus, u° G 

We make use of standard definitions for the norm and the semi- norm in H^{Q): 

\\u\\HT-(n) = {\\u\\o,n + W\m{Q)) ^ ' Hh^u) = l|Vn||o,Q. 

Since H^/'^{d^) is the trace space of H^{yi), the norm and the semi-norm in H^/'^{dVl) can be 
equivalently written as follows 

U\an=u 

Hm/2^gn) - inf PlH^n) = l|Vti||o,o- (3.6) 

U\aa=u 

Now we can define equivalent norms in H^^'^(K) and H^^'^{K): 



l^lli?i/2(ii-) — \u lHi/2(on) 



(3.7) 



U\k=u 
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where || • \\fji/2(^f^-^ is defined as in (j3.7p . because K is an open surface. 

From (j3.7p one can easily derive the expression for the corresponding H^/'^{K)-i\meT product. 
In fact, applying the parallelogram law twice, integrating by parts, and recalling notations 3°, 
4°, we find (see also 



{u,v)jj 



l/2(ft') 



o,n \ on / o,dn 



I 



du° 



\ 9n ' / o,K \ ^ dn / o,K 



(3.9) 



The space H-^/^{K) is the dual space oiH^/'^{K). We prove the following result regarding 
an equivalent norm in 

Lemma 3.1 For any u G H~^/'^{K) there holds 

||n||^_i/2(;^) ~ ||Vui^||o,n- (3.10) 
The H^^/"^ -inner product corresponding to the norm on the right-hand side of ()3.10p reads as 
{u,v)jj-i/2(^K) = {u,Vk)o,K = {uK,v)o,K Vu, G F"^/^ (i^T) . (3.11) 

Proof. Using notations 2° — 4°, we integrate by parts to obtain for any u £ H~^/'^{K) and any 



I duK 



' o,n \ dn 
Hence, we find from (j3.ip and ([3 



duK ~\ , /duK ~\ , . 

o,dn \ on / o,K \ on / o,dfi\K 



sup 



0,Q 



07^1,6^1/2 (if) ll^lli/i/2(i^) 



sup 



(3.12) 



o,n 



Let w := uk\k- One has w G H^^'^{K) because uk = on di}\K. Moreover, w° = UK\dn and, 

due to the uniqueness of the solution to the Dirichlet problem (|3.5p . we conclude that w° = uk- 
Therefore, 



sup 



> 



o,n 



iVui^lloo- 



(3.13) 
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On the other hand, it is easy to see that 



sup 



< i|Vui^||o,n. 



(3.14) 



Now ([31^ immediately follows from ([3l^ - (f3l^ . 

Using (jS.lOp together with the parallelogram law we find 



{u,v)h-i/2^K) 

Hence, integrating by parts and using notation 2°, we derive p. lip . 



□ 



The following lemma states an analogous result for the space H ^/"^{K) which is the dual 
space of H^I^{K). 



Lemma 3.2 For any u (z H ^^'^{K) there holds 

\\'^\\h-'^/2{k) ~ 



K 



(3.15) 



The H ^/"^ -inner product corresponding to the norm on the right-hand side of (|3.15p reads as 

(n,^;)^_V2(^) = = ((^)/^'^)o^ ^^'^ ^ H-'/\K). (3.16) 

Proof. Let u G H~^/'^{K). Then u° G H~^/'^{k) C H-^/'^{k). Using <^ and <^ we have 



\u 



sup 



sup 



0^w£H^/2{K) \\^\\m/'2{K) O^weH^f^iK) W^Wm/'^iK) 

\{u,V)o,k\ 



sup 



sup 



V\k=v 



sup 



\{u,v) 



0,K\ 



inf |jy||^i/2(^) 



V\k=v 



sup 



\{u,v)o,k\ 



|U,||^-l/2(^). 



Hence, using (j3.10p with u replaced by u° and with K replaced by K, we prove (j3.15p : 



FIIh-1/2(A') - \\U ll//-l/2(i^) 



Then, applying the parallelogram law, integrating by parts, and making use of notations 2°, 4°, 
we derive p.l6p . □ 
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Remark 3.1 The same arguments as above can he used to find equivalent norms and corre- 
sponding inner products in the Sobolev spaces on edges ii C dK. In particular, we will need the 
H^^'^{ii)-norm and corresponding inner product. Using the notation analogous to 3° and 4°, we 
have (cf. ^) 

Vn G H^/'^Ui), 

The next lemma states the fact that for a constant function v in (j3.16p the i?^^/^(i^)-inner 
product reduces to the L^(i^)-inner product. 

Lemma 3.3 For any u € H^^^^{K) there holds 
Proof. We have by (f3l6]) 

{u,l)jj^if.2(^K) = {u,ip\K)o,K, (3.17) 

where ip{x) (x = {xi,X2,X3) S $7 = x (0, 1)) solves the following mixed problem (see (j3.16p 
and notations 1°, 2°, 4°): find if G i7^(Q) such that 

Av9 = 0inl7, ^ = 1 on Ti = K, |^ = on (i = 2, . . . , 5), = on Tg. 

It is easy to see that cp = 1 — X3. Then iflx = 1x3=0 = 1 and the assertion follows from ()3.17p . 

□ 



3.3 Auxiliary lemmas 

The Laplace-Beltrami operator Ar will be a useful tool in our analysis. The following two 
lemmas establish its regularity separately on closed and open piecewise plane Lipschitz surfaces. 
For proofs we refer to [15| Theorem 8] and \12\ Proposition 4.11], respectively. 

Lemma 3.4 Let T be a closed Lipschitz polyhedral surface, ip € HliT) for s > —1, and let 
(f) E H^{T)/C be the unique solution to the problem 

(Vr0,Vr'^) = (^,(^) V0GFi(r)/C. 

Then cp G H^+'^iT) and 

\\4>\\m+^{r)/c < C\\ip\\H={r) 

for any r < min{s*, s + 1}, where s* > depends on the geometry ofV in neighbourhoods of its 
vertices. 
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Lemma 3.5 Let T be an open piecewise plane Lipschitz surface. Let tp S H^{r), s > —2 (resp., 
ip S -^*(r), —1 < s < (■0, 1) = 0) and (f) G H^{T)/C be the unique solution to the problem 

(Vr 0, Vr 4>) = 4>) V0 G (r)/C. 

Then (j) G H^+''{V) and 

\\(t^\\m+-r(r)/c <C\\'4'\\H-{r) (resp., ||0||//i+r(r)/c < C* llV'|[j7s(r)) 

for any r < min {s*, s + 1}, where s* > depends on the geometry ofV in neighbourhoods of all 
vertices off. 

Remark 3.2 If s + 1 < s* in the above two lemmas for given T and s, then both results are 
valid for < r < s + 1 . 

In the following lemma we formulate some useful properties of the Piola transform. 

Lemma 3.6 Let and K be two open subsets of such that = T{K), where T is 
a sufficiently smooth one-to-one mapping with a sufficiently smooth inverse 
Assume that discmK^ ~ pj^u ~ h, diamiC ~ pK — 1, and the mappings T, satisfy the same 
relations as in (|2.ip - (j2.2p . Let (p and c[ be a scalar function and a vector field, respectively, 
defined on K , and let ip = (poT^^, q = A4{q) be defined on (here, M is the Piola transform 
associated with T, see ()2.3p ). Then 

(9J,divq)o^^h = ('^,divq)o,K, (3.18) 

l|q|lo,A'h - llq||o,/f (3.19) 

if p> G L'^{K) and q G H(div, K). Moreover, for s G [0, 1], there holds 

l|q||H=(i^) < c llqllH={i^'') (3.20) 

ifqeWiK^); 

l|div q|| j^-s(^) < C h^'^ ||div q|| jj-.(^h) (3.21) 

ifdi\c\£H~'{K^); 

||divq||j^-.(^) ~ h^-' ||divq||^-s(^;,) (3.22) 

i/divq G H-'{K^). 

Proof. Relations ()3.18p . (|3.19p are well-known (see, e.g., Lemmas 1.5, 1.6 in Chapter III 
of ^). If q G il^{K^) then (cf. I30j) 

llq|lHi(i^) - llqllHi(ii'ft)- (3.23) 

Then ()3.20p follows from ()3.19p and ()3.23p by interpolation. 
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In order to prove (j3.2ip and (|3.22p we use ()3.18p and the standard scaling argument for 
scalar functions. For instance, in the former case we have 

< C sup , ^ II II = Ch \\d\Yc\\\^_s,j^hY 

The proof of (j3.22p is analogous. □ 

The following lemma states the inverse inequality for polynomials on K. We refer to [23] for 
a proof. 

Lemma 3.7 Let Vp S Vp{K). Then for any s, r ^ [—1, 1] with s < r there holds 

\\v\\Hr{K) < Cp^^''"'') MlH-iK), 

where C is a positive constant independent of p. 



4 Decompositions 

The main tool in the analysis of the EFIE is the Helmholtz decomposition of the energy space 
X. It is used to prove an inf-sup condition for the electric field integral operator and to establish 
the unique solvability of the EFIE on F (see, e.g., [151 U.2\). The following statement establishes 
the Helmholtz decomposition of X on a (closed or open) Lipschitz polyhedral surface F. This 
result has been proved in [14, Theorems 5.1 and 6.4] (for open surfaces see also [121 Section 2.4]). 

Theorem 4.1 Let 

W = {w € X; divr w = 0}, 
V = {v G X; (v,w) = Vw E WnL?(F)}. 

Then there holds 

X = VeW. (4.1) 

Furthermore, V and W are closed subspaces of X, and they can be written as 

V = Vr7^(F), W = curlr(F^/2(p)/c^ 

ifT is a closed surface, and 

V = Vr7i(F), W = curlr F^/^(F) 

ifT is an open surface. 
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In this paper we discretise the EFIE by the /ip- version of the Galerkin BEM based on the 
sequence of the RT-subspaces X/^p C X (see ()2.4p . ()2.5p ). To prove the weh-posedness of ()2.5p 
(see Theorem 12. ip we fohow \12\ 118] and consider L^(r)-orthogonal discrete decompositions of 
X/jp mimicking the Helmholtz decomposition of X: 

Xftp = Vhp e Whp, (4.2) 

where 

W/jp := {vi^ € X/,p; divr vi^ = 0}, 

Yhp := {vEX/,p; (v,w) =0 VwG W;,p}. (4.3) 

It is easy to see that W/^p C W, however, in general, V/jp ^z! V. That is why the discrete inf-sup 
condition (and thus, the unique solvabihty of (|2.5p and quasi-optimal convergence of the BEM) 
cannot be deduced by standard arguments, which are usuahy apphed to conforming Galerkin 
discretisations of coercive variational problems. 

Sufficient conditions to establish the well-posedness of the Galerkin BEM applied to problem 
p.ip were found in [12]. It turns out that it is enough to prove that discrete decompositions 
(14. 2p are in some sense close to the Helmholtz decomposition (14. ip of X when the dimension of 
the discrete space tends to infinity. The abstract formulation of this approach is given in the 
next theorem (here, we quote [18\ Theorem 4.1], see also Proposition 4.1, Corollary 4.2, and 
Theorem 4.5 in |12j). 

Theorem 4.2 Let {X„}„ be a sequence of closed subspaces X„ C X with decompositions X„ = 
V„ © W„ which are stable with respect to complex conjugation and which satisfy the following 
assumptions: 

(Al) the family {X„}„ is dense in the space X, namely 

|JXn = X; 

n 

(A2) the spaces V„ and W„ are such that W„ C W and 

sup inf ^^-^ — — >■ as n ^ oo. (4-4) 
v„ev„\{0} vev ||vn||x 

Then there exists uq such that for all f E X' and n > uq the Galerkin system 

a(u„, v) = (f, v) Vv G X„ 
has a unique solution u„ G X„ which converges quasi- optimally, i.e., 

||u - u„||x < Cinf{||u - v||x; v G X„}, 
where u G X zs the solution of (jl.ip . 
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It has been proved in [12t 118] that L^(r)-orthogonal discrete decompositions mimicking the 
Helmholtz decomposition of the space X satisfy assumptions (Al), (A2) of Theorem 14.21 with 
respect to the mesh parameter h, i.e., in the framework of the /i-version of the BEM for the EFIE. 
To prove this result for the /ip- version on quasi-uniform meshes, we will need an H~^/^(div)- 
conforming p-interpolation operator, which is introduced and analysed in the next section. 



5 Interpolation operators 

The main purpose of this section is to introduce and analyse a new H~^/^(div)-conforming 
p-interpolation operator. This operator is necessary to deal with low regular vector fields, 
such as gradients of solutions to boundary value problems for the Laplace-Beltrami operator 
on polyhedral surfaces (see the regularity results of Lemmas 13.41 and 13. Sp . We will construct 
the H~^/^(div)-conforming interpolation operator by employing the ^"^/^-projection for the 
divergence term. For this operator we then prove quasi-stability with respect to polynomial 
degrees and commutativity with the ^^^/^-projector. 

In this section we use standard differential operators div, curl and V, curl acting on 2D 
vector fields and scalar functions, respectively. First, let us recall some known interpolation 
operators acting on scalar functions and vector fields on K. Let 11^ : L'^{K) —i- Vp{K) be the 
standard L^-projection onto the set of polynomials Vp{K). We will also use the i7^^/^-projector 
onto Vp{K) denoted by Hp^^^ : H-^/'^{K) Vp{K) and satisfying 

(u-n-i/2n,?j)_i ^ = yv£Vp{K). 



2 ' 



Here and below (•,•)„! ^ denotes the H ^/^(i^)-inner product (see (j3.16p ). 

In [21] two projection-based interpolation operators have been introduced and analysed. 
These are the if ^-conforming interpolation operator Up : H^^'^{K) Vp{K) and the H(curl)- 
conforming interpolation operator 11^'^''^ : 't¥\K) fl H(curl, ii') 7-'p'^'^{K) (here, r > in both 
cases and Vp'^'^{K) = Vp-i^p{K) x Vp^p-i{K) is the H(curl)-conforming (first) Nedelec space 
of degree p). Later, in [S] these operators were employed to prove the discrete compactness 
property for hp adaptive rectangular edge finite elements. Due to the isomorphism of the curl 
and the div operator in 2D (and, as a consequence, the isomorphism of the Nedelec elements of 
the first type and the RT elements), one can use the results of [TP, ^ related to the operator 
IIp'^''^ in the H(div)-settings. 

We will denote by IIp'^''' the corresponding H(div)-conforming projection-based interpolation 
operator. Then for r > the following diagram commutes (see Proposition 3 in |21j): 

H^+r^K) ^ H'^(i^)nH(div,i^) ^ L\K) 



n°_i (5.1) 



VpiK) ^ Vf'iK) ^ Vp.,{K). 

Furthermore, we use the above mentioned isomorphisms to reformulate the following two results 
from p]. 
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Lemma 5.1 [H Theorem 4] Let Ap = Vf^'^{K) and Bp = curlP°(K) © V div7']^T'°(K). 

Then the following stability condition holds 



aeApbeBp ||a||o,i^ ||b||o,i^ \{p + l){p + 2) 

This lemma and the definition of the interpolation operator IIp"''^ imply the following L^- 
stability result for the version. 

Lemma 5.2 [9, Theorem 7] Let u € H''(i^)nHo(div, K), r > 0, be a curl-free bubble function on 
K, and letuf'' := Up"'°u. Thenuf^ G Vf^'^{K) is discrete curl-free (i.e., {uf^ , curie/)) q^k = 
for any cj) € Vp{K)) and there holds 

||u - Up"'||o,_ft' < Cp^/^ inf ||u - qp||o,_ft:. 

In the next lemma we estimate the error of the best L^-approximation of low regular vector 
bubble functions by RT-elements of degree p. 

Lemma 5.3 Let u E if^K) f] Ho(div, A'), r € [0, 1], be a general bubble function on K. Then 

inf qpWo^K < C p"'' \\u\\iirrKy 

Proof. For r = 1, the statement follows from |9j Lemma 8]. For r = it is trivial, because 

inf ||u - (ipWo K < llullo,^. 

Then we obtain the whole range r G [0, 1] via interpolation. □ 

We will use the above two lemmas to prove the following auxiliary result. 

Lemma 5.4 Let Vp G 7-*^'^'°(i^) be such that (vp,Wp)o,ii- = for any Wp G Wp{K) := {w G 
Vf^'°{K); divw = 0}. Then there exists v G VH{K) such that 

llv||o,_ft: < C ||divvp||^_i/2(^) (5.2) 

and 

l|v - Vp||H(div,i<') < C'P^" l|divVp|lj^_i/2^^) (5.3) 
for any eo > with C = C(eo) > 0. 
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Proof. For given Vp, we solve the Neumann problem to find / G H^{K)/C such that 

(V/,V5W = -(divVp,9)o,i^ ^geH\K)/C. (5.4) 

Then we set v := V / . One has 

divv = A/ = divvp G (5.5) 

Hence, / € UiK), v € Vn{K), and ^ holds, because A : H'^{K)/C H-^{K) is an 
isomorphism. 

Note that divvp E ^^^^^(i^) for any e > 0. Therefore, the standard regularity result for 
problem (j5.4p reads as (see, e.g., [22j and cf. Lemma 13.51 and Remark 13. 2p : / E H^~^'^{K) and 

\\f\\H^+^{K)/c < l|divVp||^_i/2+E(i^) 

for any < r < ^ + e and for arbitrary e € (0, 1]. Then, using the continuity of the gradient as 
a mapping H^~^'^{K) — > tV{K), we have 

\Mn-iK) <C\\f\\m+-(K)/c<C\\divVp\\j^-i/2+e(^K), 0<r<l/2 + e, eG(0,l]. (5.6) 

Since v € H^{K) Pi Ho(div, K), we can apply the interpolation operator IIp'^''^ to define v^^^ := 
Hp'^''^ V € V^^'^iK). Recalling that 11^'^'° commutes with the L^-projector (see (|5.ip ) and using 
()5.5p . we find that 

divVp'^ = divvp = divv. 
Hence, (vp — "v^") G Wp(i^). This fact implies the relations 

(v, Vp - ^r'^")o,K = (V /, Vp - v^*^)o,i^ = and (vp, Vp - v^*^)o,i^ = 0, 

where the latter equation holds by assumptions on Vp. Therefore, 

||v-Vp||o,i^ < ||v-v^-||o,/r. (5.7) 

Since v is curl-free, we apply Lemmas 15.21 15.31 and then use inequality ()5.6p . As a result, we 
obtain 

||v-vf||o,i, < inf ||v-qp||o,i, 

< Cpi/2-'-||v||H.(i^) <Cp^/'-n|divVp||^_V.+.(^) (5.8) 
for any e G (0, 1]. Then making use of the inverse inequality (see Lemma l3.7p we estimate 

l|divVp||j:^_i/2+.(;^) < C/^||divVp||^-i/2(^) < Cp^^||div Vp||^_i/2(;^). (5.9) 
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Now we again use ()5.5p and then put together ()5.7p - (|5.9p . We obtain 

l|v - Vp||H(div,K) = ||v-Vp||o,i^ < C//^+^^"''||divVp||^_i/2(^). (5.10) 

Given an arbitrary eq > 0, we select e = minjeo, !}• Then (j5.3p follows from (|5.10p by taking 
r = i + e. □ 

H^^/^(div)-conforming interpolation operator. Now we proceed to the main goal of 
this section. Given a vector field u G H''(i^) n H^"'^/^(div, with r > 0, we construct an 

interpolant u^* = Hp"' € 7'p^(i^). In particular, u^' is defined as the sum of three terms: 

uP = ui+uf + uf. (5.11) 

The definition of ui and U2 follows the construction from [21j . Let ui be a lowest order inter- 
polant defined as 

4 „ 

ni = Y.U u-nda)cl>i, (5.12) 

where n denotes the outward normal unit vector to K, and 0j (i = 1, . . . , 4) are the standard 
basis functions for VY^{K), defined by 



0i • n 

For any edge £i C dK one has 



1 on ^j, 
on dK\e^ 



J{u-ui) -nda = 0. (5.13) 

Hence, there exists a function ip, defined on the boundary dK, such that 

= (u — ui) • n, "0 = at all vertices. (5.14) 
Then we define if) 2 G "Ppi^i) by projection 



0/ 



(V^|,,-V^<A)/^i/2(,,) = o ycpev'^ii.) (5.15) 

(see Remark 13. II for the expression of (•, •) /^i/2 (£.-)). Extending tp2 by zero from £{ onto dK (and 

keeping its notation), we denote by ■ip2p ^ ^p{^) ^ polynomial extension of 'ip2 from dK onto 
K, i.e., 

4]^erp{K), 4:pk = 4'^ 4:p\dKv, = o. (5.16) 

Then we set 

4 

"2 = X]^2A' where u|_^^ = curl ^Jp. (5.17) 

i=l 
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The interior interpolant Ug is a vector bubble function living in 7^^'^'^{K) and satisfying the 
following system of equations: 

(div(u - (ui + + )), divv)_i/2,ii- 

= (div(u-(ui+u^)),divv)_i/2,i,=0 Vv G P^T,0(^)^ (5^18^ 

(u - (ui + + u|), curl (P)o,K = V(/. G V^iK). (5.19) 

We note that any polynomial extension satisfying (|5.16|) can be used for the construction of the 
edge interpolant Ug. Nevertheless, the interpolant = ui + U2 + Ug is uniquely defined, which 

follows from (j5.17p - (|5.19p . It is also easy to see that Hp"' ^ preserves polynomial vector fields, 

i.e., Hp"' "2 Vp = Vp for any Vp G Vf^iK). 



Remark 5.1 In contrast to the Lp'-inner product employed to define the H^cml)- conforming 
interpolation operator (and its 'H.{div)- conforming counterpart) in [21j, we use the H~^^^-inner 
product in (|5.18|) . 

Proposition 5.1 For r > the operator 

uf^'~^ : W{K) n ii-^/'^{dw, K) ^ i?{K) n ir^f'^{d\Y, k) 



is well defined and hounded. For arbitrarily small £9 > the norm of this operator satisfies 

< Cp''>, 

c 

where C > is independent of p but depends on Eq and r, and \\ ■ \\c denotes the operator norm 
in the space c{w {K) r\il-^''^{di-v, K),\?{K) r\n-^/'^{di-v ,K)^ . 

Proof. Let u G (ivT) nH~^/^ (div, K)^ r > 0. We will study each term on the right-hand side 
of (jS.lip . Throughout the proof we denote by s a small parameter such that < s < min r} 
for given r > 0. 

Step 1. Fixing an edge C dK and using a function 



1 on ^j, 
on dK\ii 



as a test function, we integrate by parts to obtain 

J u ■ nda = J {u ■ n) (pi da = J (div u) (pidx + J u • V0j dx 



dK K K 

< ||divu||jj_i+3(j^) \\(l)i\\m-s{K) + l|u||H=(i^) IIV^illH-'-Cif) 

< C{(l)i,s) niullHr(ii') + ||divu||^_i/2(^) 
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Note that if divu € H^^^^{K) then an extension to divu € H^^^^{K) exists but is not unique. 
However, by assumption divu € H~^^'^{K) C H^^^'^{K), which is a unique extension (see [27j 
for details). Thus, ui in (j5.12p is well defined. Moreover, since ui is a lowest order interpolant, 
we find by the equivalence of norms in finite-dimensional spaces that 

lluillHCdiv.A-) ^ C'X] I / " ■ "^'^ - ^ (|I"IIh'-(e-) + l|divu||^_i/2(^)^ . (5.20) 

Let use denote by 7^^^"*^ a right inverse of 7tr with 7^""^ : H^^'^~^ (dK) H^^'^(K). Taking an 
arbitrary v € H^^^~^{dK) we integrate by parts similarly as above to estimate 

y (u — ui) ■ n V da 
dK 

< ||div(u - ui)||j^_i+,(^) \htT^v\\Hi-s(^K) + ||u - uiIIh.(a') II V(7tT^^)llH-(;^) 

< C (^||u - uiIIhs(^) + |[div(u - ui)||^_i+,(^)^ ||f ||j:^i/2-s(9^). 

Hence, (u — ui) • n G H^^/^^^{dK) and, due to the finite dimensionality of ui, we obtain by 
using estimate (j5.20p : 

ll(u-ui) -nll^-i/a+s/gj^) = sup — 

< C(^||u-ui||h»(;^) + ||div(u- ui)||j^_i+,(^)^ 

< C (^\\u\\u=(^K) + l|divu||^_i+,^^^ + ||ui||H(div,/^)^ 

< C (^llullH'-(ii-) + ||divu||^_i/2(^)) . (5.21) 
Step 2. From the construction of ui and from the result of Step 1 we conclude that 

(u - ui) • n G H-^l^+\dK), J (u-ui) -nda = 0. 

dK 

Therefore, due to the isomorphism ^ : H'^/'^+'{dK)/C H^^^'^'^' (dK) (see j[21i Lemma 2]), 
the function ip in (f5T4]) is weh defined, ip G H^/'^+'idK), ip\i,. G H'^^^iii) for any edge ii C dK, 



and 
4 



5ZllV'l4lli?i/2(4) ^ 2^ II^I^J/fi/2+^{4) <CM\Hi/^+s(dK) < C \\{U - Ui) ■ n\\fj-i/2+s(^9K)- 

(5.22) 



i=l 1=1 
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Hence, (|5.15p is uniquely solvable and 

\\Tp2 11^1/2(4) < C lli^l/2(4)- (5.23) 

Furthermore, applying the polynomial extension result of Babuska and Suri j^, we find the 
desired polynomial ^ T^pi^) (see (|5.f6p ) satisfying 

\\4:p\\miK)<CU2\\m/^ie,y (5-24) 
Thus, in (ICTT]) is weh defined. Putting together we find 

4 4 
||u^||o,A- < ||curl^Jp||o,i^ < |[V'Jp||//i(K) 

i=l 1=1 
4 4 

< <^5^ IIV'2 lli^i/2(£,) < llV'kJli/i/2(4) < C||(u-ui) •n||^_i/2+.(g;^). 

i=l 1=1 

Hence, making use of (|5.2ip . we obtain 

l|u2llH(div,K) = lluallo.ii- < C (||u||Hr(A-) + ||divu||j^_i/2(^)) . (5.25) 

Step 3. The vector bubble function Ug is uniquely defined by (j5.18p - (|5.19p . To estimate 
the norms of and divu^ we use the discrete Helmholtz decomposition ()4.2p restricted to 
X,{K) = '°(i^) 

u| = Vp + curl(/)p, (5.26) 

where (/>p G Vp{K) and Vp G 7-*p (i^) is such that {'Vp,Wp)o^K = for all Wp G Wp{K) (see 
Lemma 15.41 for the definition of Wp{K)). 

From (|5.18p one has by using the result of Step 1 



||divu^||j:^_i/2(^) < C ||div(u - ui)||^_i/2(^) < C (^||div u||j:^_i/2(^) + |divuij 

< C (^llullH'-(ii-) + ||divu||^_i/2(^)) . (5.27) 

Then we apply Lemma 15.41 there exists v G V7i{K) and an arbitrarily small eo > such that 

||vp||o,/r < ||v - Vp||o,i^ + |Iv||o,_ft: < Cp^o ||divvp||j^_i/2(^). 

Hence, recalling that divvp = divug, we obtain 

hp\WK < Cp'^' ||divuP||^_i/2(^). (5.28) 

Since (vp, curl(^p)o,_ft' = 0, we estimate the norm of curl (j)p by using (j5.19p and by employing 
the results of the first two steps: 

||curl0p||o,E- < llu - ui - u^\\o,K < C (^llullH'-(ii-) + ||divu||^_i/2(^)^ . (5.29) 
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Combining (|5.27p - (|5.29p and applying the triangle inequality we obtain by making use of de- 
composition (|5.26p 

llusllo.A' + l|divu?^||j^-i/2(^) < Cp^" (||u||H'-(i^) + ||divu||^_i/2(^) 

To finish the proof it remains to combine the results of the three individual steps and to apply 
the triangle inequality to decomposition (jS.lip . □ 

Proposition 15.11 proves a quasi-stability of the interpolation operator Hp"' ^ . The following 
proposition states the commuting diagram property for IIp'^' ^. 

Proposition 5.2 For any u G W{K) f] H^-'^/^(div, K), r > 0, there holds 

div = Hp^f (divu). (5.30) 

Proof. For any (f € Vp-i{K) there exists Vp G 7^^'^{K) such that divvp = (f. Therefore, 
decomposing IIp'^' as in (j5.1ip . we need to show that for all Vp € T^^'^^K) there holds 

div^u — IIp'^' ^u^,divvp^ ^ = (div(u — (ui + Ug)), divvp)_x/2,_ft' = 0- (5.31) 

.div,-i 



Let us also decompose Vp = Hp"' ^Vp G 7^^'^{K) as in (|5.1ip : 

Vp = vi + + vf , div vi = const, div = 0, vf G Vp^'°{K). 
Then, recalling (j5.18p . applying Lemma 13.31 and integrating by parts, we prove ()5.3ip : 



div ( u — IIp'^' ^ u ) , div v 



1/2,K 

(div(u — ui — Ug), const)_i/2,A' + (div(u — IIp'^' J , divvj 



-1/2,K 



= (div(u — ui — Ug), const)o,i<' = const J (u — ui — Ug) ■ nda = 0. 

dK 

For the last step we used the fact that Ug • n|gx = and then applied (j5.13p . □ 



6 Proof of Theorem 12.11 

In this section we prove Theorem 12.11 relying on the abstract convergence result of Theorem 14.21 
Since the discrete decomposition (|4.2p is stable with respect to complex conjugation, one needs 
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to check that assumptions (Al) and (A2) are satisfied. First, we note that the family {X.hp} of 
RT-spaces is dense in X*^. Since the injection X° C X is dense as weU (see, e.g., [171 Lemma 2.4]), 
we conclude that the family {X/jp} satisfies assumption (Al) of Theorem 14.21 

It was mentioned in Section U] that W/ipCW by construction. Thus, it remains to prove that 
the subspace Yhp defined by (|4.3p satisfies assumption (j4.4|) . In particular, we will show below 
that there exists a sequence {5hp}, 6hp 0, such that for any given v^p S V/jp there exists 
V E V satisfying 

||vhp - v||x < (J/ip||v/ip||x- (6.1) 

We prove ()6.ip for a closed (resp., open) surface T. The proof consists of 5 steps. 

Step 1: Construction of v. Given v^p € V/^p, we solve the following problem to find 
/ G H\r)/C such that 

(Vr/,Vr5) = -(divrv;,p,g) G H\r)/C. (6.2) 

We set V := Vr /. Then 

divr V = divr Vhp G L'^{T), (6.3) 

and, due to Theorem 14. H there holds v G V. Note that divpv/jp G H~^^'^{T) (resp., divp v/jp G 
H~^/'^(r)) and (divp v/^p, 1) = 0. Therefore, the regularity result for problem (|6.2p reads as (cf. 
Lemma [3^ (resp.. Lemma [33]) ): / G H^'^^{T) with r = min{s*, ^} — ei, where s* > depends 
on the geometry of F, ei > is arbitrarily small. Moreover, using the continuity of Vr, we 
conclude that v G H^_(F) and 

l|v||H'-_(r) < C ||/||Hi+'-{r)/c < C ||divr v/jp||^_i/2(r), r = min{s*, ^} - ei (6.4) 

(here and below we use the convention H~^^'^{T) = H^^^'^(r) if F is closed). 

Step 2: Bounding ||v/ip — v||x by ||v — v/ip||o,r' In view of (j6.3p we can estimate the norm 
on the left-hand side of (16.11) as 



hhp - V||x = II Vhp - V||jj_i/2^p^ < C \\vhp - v||jj_i/2+.i^^^ 

= G sup J — II (6.5) 

weH^Z^-=i{r)\{o} l|w|lH^/2-=i(r) 

with the same ei > as in (fO]) . Let w G H^Z^"^'(F). Then, by [H Theorem 3.4] (resp., O 
Theorem 6.1]), there exists a unique pair , W2 G H^{r)/C (resp., wi G -ff"'^(F)/C, W2 G i7o(F)) 
such that 

w = Vr t«i + curlr ^2 • (6.6) 
Moreover, W2 is the solution of the problem 

— Ar ^2 = curlr curlr W2 = curlr w. 
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Since curlpw G H^'^^'^^^'^iT) and (curlpw, 1) = (resp., curlp w € i/^^/^^^^(r)), we apply 
Lemma 13.41 (resp., the Dirichlet analog of Lemma l3.5p to prove that W2 € //^^'(F) (with the 
same r as in (16.41)1 and there holds 



\w2\\m+r{T)/c < llcurlr w||j^-i/2-.i(r) ^ C \\w\\^i/2-ei^^^ 
(resp., \\w2\\m+r(j^) < C ||w||jji/2-.i^j^p. 



(6.7) 



Therefore, W2 G H^^'^iJ'j), r > 0, for any element Tj, and one can apply the if ^-conforming 
interpolation operator 11^ to find a continuous piecewise polynomial w^^ defined on F such that 
^2 j ^2^\^j ° = ^pW2,j £ 'Pp{K) (here W2j := W2\tj ° Tj). To estimate the iJ^-semi-norm 
of the error {w2 — we apply the standard scaling argument (cf. [191 Theorem 4.3.2]) and 

the approximation result for Hp on the reference element (see \21[ Theorem 2]): 



\w2-W2^\Hl{rj) < C 



c 



W2,j - ^lw2,j 



m(K) 



ipp - It{w2,j - ipp] 



HHK) 



< Cp ^ inf \\w2,j - ^p\\m+r(K), < £2 < r < 1/2. 

Let s = 1,2. Using Theorem 3.1.1 of [19] and the scaling argument, one has 



(6.8) 



Therefore, by interpolation, 



inf 



\W2,j - 'P\\h''(K) 



< C\w2,j\H'>iK) < C h" ^ \\W2\\H''{TJ)■ 



- ij^f , J I ^2, i - 1 1 1+- (ii-) < C' /l'' 1 1 ''«2 1 1 !+'■ (r, ) , 



and from 



we conclude that 



\w2-w!^^\Hi(r,)<Ch''p \\w2\\H^+r(^^J)■ 



Ilence, for a closed (resp., open) surface F we obtain 

1/2 



I hp I 

\W2 - W^'^Ih^T) 



Ek2-^«2''lii(r,) 



\w2\\H^+r(r)/c 



(6.9) 



(resp., \w2 - W2''\H^r) < C p ^2) ||^t;2||_f/i+'-(r))- 



In addition, for the case of an open surface, we note that vanishes on dV. Then recalling 
the commuting diagram property for 11^ (see ()5.ip ) and the definition of X/^p (see ()2.4p ). we 
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conclude that curlru;2^ € X/^p C X. Moreover, curlp 1x^2^ £ W/^p C (W fl L^(r)). This fact 
together with the L^(r)-orthogonahties V ± (W n L^(r)) and V^p 1- ^hp imphes 



(v - Vftp, curlr w!^'^) = 



(6.10) 



Now we use (|6.3p . (j6.6p . (|6.7p . (j6.9p . (j6.10p (and (j3.4p . if F is an open surface) to prove for any 



w G H 



l/S-ei 



(r) 



(v-v,,p,w) = (v - v,,p,Vr'u;i) + (v - v/jp, curlr ^112) 

= -(divr(v - ^^hp),wi) + (v - v/jp, curlr (^2 - ^^'2'')) 



< V 



hp I 



v/ip||o,r |^^2 - ^^2 li/i(r) 



nrn < C/l"p-("-^2) llv 



v/ip||o,r ||w|| 1/2-ei 



^(r)- 



(6.11) 



P 

Using this estimate in (j6.5p we find 

||v/ip - v||x < C/i''p"(''"^2) ||v _ v/ip||o,r, 

where r = min{s*, ^} — ei, ei > 0, £2 £ (0;'^) (^ii ^2 can be arbitrarily small). 

Step 3: Bounding ||v - v/jp||o,r by ||v - uf^' ^^V||o,r. We recah that v G H'^_(r) n 
H(divr,r), r S (0, ^). Therefore, vlr^ € H'^{Tj) n H(divr,rj) for each element Tj and we can 
define v^^ € X/jp such that 

div,-i , 



A^-^(v|r,) , 



where A4j is the Piola transform (see (j2.3p ). 

Using commutativity (|5.3Up and the fact that divr v = divr v/^p is a piecewise polynomial 
on r, we find 



divr v^p = divr v = divr v?,p. 



Hence, (v^p-v^^ ) G W/^p C (WnL^(r)), and using again the orthogonalities V ± (WnL|(r)), 
V/jp ± Whp we have 



(v - ^hp,^hp - v; 



div\ 
hp I 



0. 



Therefore, 



|v - v/,p||o,r < II V - ||o,r- 



(6.12) 



_div, 



Step 4: Estimating ||v — v^JjUo^r- Using (j3.19p and the quasi-stability of Hp ' ^ (see 
Proposition 15. ip . we estimate 



|Vhp |r,||o,rj 



MAlif^-HMj\^\T^)) 



o,r, 



< C 



0,K 



< (||A^7^(v|r,)llH'-(i^) + ||div(>lTi(v|r,))||^_v2(^) 

< (||X7^(v|r,)||H'-(i^) + l|div(X-i(v|r,))||^-i/2+.3(i^)), 
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where £0)^3 > are arbitrarily small. Hence, recalling that div(A4j ""^(vlrj)) = div(A^j ^{'^hplvj)) 
is a polynomial on K, we make use of the inverse inequality (see Lemma l3.7p and then apply 
(I32QD, (13:211) : 

||v^^-|r,|lo,r, < Cp^o+2e3 (^\\Mj\v\r^)\\nriK) + \\diviMj\v,,\rMH-yHK) 
< Cp''+^'' (||v||H'-{r,) + /i'/' \\divr^hp\\H-y^ir 
Therefore, 

/ \ 1/2 

||v-v^-||o,r < Mo,r + i^\Wf;\rX,T,) 

j 

< l|v||o,r + C/»+2^« rj;(||v||2j.(p^,) + /i||divrv;,p||^_i/2(r^.))) ^ 

i 

< Cp^o+2e3 (|[v||H._(r) + lldivr v;,p||^-i/.(r)) . (6.13) 



Step 5: Conclusion. First, we select all small parameters {k = 0,1,2,3) such that 
£o + £i + £2 + 2^3 = e for any given e £ (0, min {s*, ^}). Then putting together the results of all 
previous steps, i.e., estimates ()6.4p and (|6.1ip - (j6.13p . we obtain 

/ \min{s*,i}-e / \min{s*,i}-e 

||vhp - v||x < C j ||divr v/jp||j^_i/2(r) < (^pj l|v/ip||x- 

This proves (16. ip . Therefore, the subspace V/ip satisfies (14. 4p . Thus we have shown that the 
discrete decomposition ()4.2p verifies assumption (A2) of Theorem 14.21 in the framework of the 
hp-version of the BEM with quasi-uniform meshes, and the proof is finished. 



7 Concluding remarks 

The main result of this paper - the convergence of the /ip-BEM with quasi-uniform meshes (and 
thus the convergence of the h- and the p- versions as particular cases) for the EFIE - is proved 
for a sequence of Raviart-Thomas spaces on quadrilateral elements, which can be parallelograms, 
convex quadrilaterals or curvilinear ones. To that end it was essential to introduce and analyse 

a new H~^/^(div)-conforming interpolation operator Hp ' ^ on the reference element. To show 

the stability of IIp'^' ^ (with respect to polynomial degrees) we relied, in particular, on the 
discrete inf-sup condition (see Lemma IS.ip . which was established in [9]. 

The case of triangular elements. We note that Lemma 15.11 is not available for the reference 
triangle. However, the corresponding result has been conjectured and numerically evidenced 
in [To] for H(curl)-conforming Nedelec elements of the second type. In our H(div)-conforming 
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settings it means that the conjectured result of Lemma l5.ll is numerically justified for Brezzi- 
Douglas-Marini spaces on the reference triangle (this is due to the isomorphism of the curl and 

the div operator in 2D). Since all remaining arguments in the construction of IIp'^' ^ and in 
the proofs of Propositions 15.11 15.21 and Theorem 12.11 are valid for BDM-spaces on triangles (cf. 
also [21j), we conclude that the main result of the paper holds for these spaces, provided that 
the discrete inf-sup condition discussed above is true. 

Regularity results and error analysis. In this paper we do not present the regularity result 
for the problem under consideration, neither perform an a priori error analysis for the /ip-BEM. 
However, we note that in [8j we derived explicit formulas for typical singularities inherent to 
the solution of the EFIE on piecewise smooth (open or closed) Lipschitz surfaces. The p- 
approximation analysis of these singularities (including the least regular ones) was performed 
in [8] on a plane open surface. This analysis relied on our results for the Laplacian, see [UlTllS], 
by using continuity properties of the surface (vector) curl operator. Since those results and 
properties (in corresponding spaces) are valid for polyhedral and piecewise plane open surfaces, 
the proof of an a priori error estimate for the p-BEM from [8] carries over without essential 
modifications to the more general case considered in the present paper. We stress that this 
p- approximation result holds only for affine families of meshes. The error analysis of the Zip- 
version on quasi- uniform meshes presents more difficulties. In particular, the involved Sobolev 
spaces of negative order are not scalable under affine transformations. Therefore, this analysis 
is not a trivial extension of our p- approximation results. 

Non-affine quadrilateral meshes. Although the main convergence result of the paper holds for 
non-affine quadrilateral meshes, it is not clear if the RT-spaces in this case would provide optimal 
approximations in the energy norm for the EFIE. For instance, in the /i-version of the finite 
element method the degradation of convergence rates (in the H(div)-norm) was observed for 
RT-elements on general convex quadrilaterals (see [1]). An alternative family of finite elements 
on non-affine meshes was introduced in [1] and was shown to provide optimal /i-approximation 
order in H(div). These Arnold-Boffi-Falk (ABE) elements can also be used in the BEM for 
the EFIE. However, the solvability, convergence, and a priori error estimation of the /i-, p-, or 
/ip-BEM with ABF-elements are open problems. 
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